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Abstract 

We demonstrate a variety of ordered patterns, including hexagonal struc- 
tures and chains, formed by colloidal particles (droplets) at the free surface 
of a nematic liquid crystal (LC). The surface placement introduces a new 
type of particle interaction as compared to particles entirely in the LC bulk. 
Namely, director deformations caused by the particles lead to distortions of 
the interface and thus to capillary attraction. The elastic-capillary coupling is 
strong enough to remain relevant even at the micron scale when its buoyancy- 
capillary counterpart becomes irrelevant. 
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2D organization of nanometer- and micrometer-sized colloidal particles at fluid inter- 
faces is a fascinating phenomenon of both fundamental and applied interest. The nature 
of micron-scale interparticle forces, especially of the attractive nature, remains a subject 
of an ongoing debate [1-9]. In most cases, the interparticle forces are isotropic and cause 
hexagonal ordering. Anisotropic interactions can be achieved when the particles are non- 
spherical [4,5]. Apparently, the surfaces of anisotropic fluids, i.e., liquid crystals (LCs), 
can also support anisotropic interactions and thus a richer variety of ordered patterns as 
compared to isotropic fluids, but the current knowledge in the field is rather limited. It is 
known that small particles might form chains decorating the surface director field [10,11] 
and that the surface air bubbles can be accompanied by point defects [12]. The behavior of 
particles entirely in the LC bulk is studied much better [13-21]. Because of the anisotropy 
of molecular interactions at the particle surface (the phenomenon of anchoring [22]), the 
embedded particle causes director distortions. The most frequently met distortions are of 
dipole symmetry; they lead to chaining of droplets [13]. Recently, an unexpected and so 
far unexplained hexagonal structures have been discovered for an array of glycerol droplets 
assumed to be in the nematic bulk [18-20] and for droplets in smectic membranes [21]. 

In this work, we use confocal microscopy to demonstrate that placement of colloidal 
particles (glycerol droplets of radius R = 1 — 10 fim) at the LC surface leads to attractive 
interactions and ordered patterns of hexagonal and chain type that depend on the thickness 
of the LC film. The attractive interactions can be explained by the elastic-capillary coupling, 
as the particle-induced director deformations distort the LC surface. 

The glycerol droplets are obtained as in Ref. [18]. A Petri dish containing a layer of 
glycerol and the nematic LC pentylcyanobiphenyl (5CB, EM Industries) on top of it, is 
kept at 50°C for 10 min, to facilitate diffusion of the glycerol molecules into the 5CB layer 
which is in the isotropic phase above ~ 35°C. When the sample is cooled down, solubility 
decreases and one observes appearance and growth of glycerol droplets. The technique 
produces droplets of a constant radius R determined by the cooling rate and the number of 
thermal cyclings. The thickness of the LC film is controlled in the range h — (3 — 100) \xm. 

2 



The 5CB film is in the hybrid aligned state: the director n is parallel to the LC-glycerol 
interface at the bottom and perpendicular to the air-LC interface at the top. 

In order to determine the location of droplets, we use the fluorescence confocal polarizing 
microscopy (FCPM) [23]. Two different dyes with separated absorption and fluorescence 
bands, fluorescein and Nile red (Aldrich), were added (0.01 wt %) to tag glycerol and LC, 
respectively [23]; they did not change the appearance of patterns. The FCPM textures of 
vertical cross-sections, Fig.la-c, unambiguously demonstrate that the glycerol droplets are 
trapped at the LC-air interface. For a better clarity, the images in Fig.la,b are taken for a 
large drop; smaller drops are also located at the interface, Fig.lc. 

A 4-6 hour relaxation at room temperature results in 2D ordered structures of droplets. 
The order is hexagonal when the LC layer is thick, h » R, h = (20 — 100) /j,m, Fig.2a,b. 
The droplet center-to-center separation r eq increases with the droplets diameter 2R, Fig. 2d; 
each data point is an average of the data collected by video-recording the patterns every 6 s 
for 10 min. The angles between the directions to the nearest neighbours changes in the range 
60° ± 8° for well-equilibrated structures. At this stage, it is hard to quantify the lattice 
symmetry more precisely as it might be influenced by slowly changing director pattern in the 
bulk. In thin films, h/R m 2 5, 3fim < h < lOfim, Fig. 3a, the pattern changes completely, 
as the droplets form chains oriented along the average horizontal projection of n. The 
length of chains (number of droplets in a chain) increases when h decreases. Importantly, 
both hexagonal and chain patterns disorganize when 5CB is heated into the isotropic phase. 

The very existence of 2D hexagonal pattern in thick films does not discriminate between 
attractive and repulsive interactions among the particles, as in the confined systems such 
as a Petri dish the ordering can be induced by purely repulsive forces. The presence of 
repulsive force is manifested clearly, as the droplets are well separated, r eq R3 (3 — 4) R, 
Fig. 2d. The following three experiments show that the interaction has also an attractive 
component. (A) We start with the Petri dish in which the 2D lattice occupies the whole 
area and then use a filter paper to wipe about 1/2-1/3 of the droplets. The remaining 
droplets restore the 2D lattice within the smaller area, with the lattice parameter close to 
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the original one. (B) We change the direction of the meniscus tilt near the Petri dish walls 
by coating the walls with a hydrophobic material (stearin); the types of the 2D ordering 
remain the same. (C) After the 2D lattice is formed in a regular manner, we add (through a 
microsyringe) a small amount (0.01% by weight) of a surfactant (cetylpiridinium chloride) to 
the nematic layer. The interparticle distances decrease; the shrinked clusters are separated 
by droplet-free areas, Fig.2c. 

What is the mechanism of the droplet attraction? The gravity-mediated capillary effect 
and the van der Waals forces appear to be insignificant, as the corresponding pair potentials 
[24,22] 

U g (r) « ^-9 2 p 2 LC ln(y); U vdW (r) « - AR — ; (1) 

&lca A 12, [r — Iti) 

are of the order of thermal energy ksT pa 4 x 10~ 21 J or less, when r is a few microns. 
Here glca ~ 3.8 * 10~ 2 J/m 2 [25] is the LC-air surface tension coefficient, pic ~ 10 3 kg/m 3 
is the LC density, g is acceleration due to the gravity, A = \J<Jlca/ (qPlc ~ ff) is the 
capillary length; II pa fT//i 2 is the elastic disjoining pressure for the hybrid aligned film [22]; 
K ~ 10 _11 ./V is the average nematic elastic constant; in all cases of interest to us, A >> r eq) 
as A ~ 2mm when h — ► oo and A ~ 0.5mm when h = 5pm; finally, A is the Hamaker 
constant of the order of 10" 19 J or smaller [22]. 

Attraction of droplets at the nematic surface (and absence of attraction and ordering 
when 5CB is in the isotropic state) can be explained by the coupling of orientational elasticity 
and capillarity. The elastic distortions caused by the submerged part of the particle give 
rise to a vertically resolved force j e \ that shifts the LC-air interface upward or downward, 
depending on the particle properties, to reduce the elastic energy of the LC host; f e i is 
balanced by surface tension. Non-flat interface causes long-range attraction of particles. 
As the first approximation, we divide the total interaction potential U (r) into two parts, 
U e ic (r) describing the elastic-capillary coupling and U e ib (r) describing the "pure bulk" elastic 
interactions. 

To estimate the orders of magnitude, we use the idea that f e i is related to the 



non- vanishing surface anchoring of n(r) at the LC-particle interface: f e i = -j^E e i(z) ~ 
§ n ^W(h ■ 9) 2 ds, where E e i(z) is the total elastic energy of distortions caused by the 
particle, Q(z) is part of the droplet surface immersed in the LC, W is the (polar) anchoring 
coefficient of the LC-glycerol interface, z> is the normal to the interface, and z is the vertical 
coordinate of some fixed point, say, the bottom of the droplet. For f ei AnWR, with 
W = 10~ 5 J/m 2 [22], one finds f el ~ l(T 10 iV ~ 10 K when R = lfim. The force f el replaces 
the buoyancy force in the standard treatment [3,24]; it is balanced by the surface tension, 
hence, f e i = 2n<jLCAri simp, where ip is the meniscus slope at the triple contact line of radius 
r\. The vertical displacement £ (r) —r\ sin ^ In ^ changes the interfacial area and leads to 
the attractive potential written for two identical particles as U e i c = — / e ;£, or 

Ueic(r) = -G^H r -), (2) 

similar to the buoyancy case [3,24]. For f e i ~ 10i^, one estimates U e i c ~ 3 x 10 _19 J ~ 80 
ksT. The elastic-capillary coupling is strong enough to keep droplets attracted even at 
micron scales where the gravity mechanism vanishes: the dependence U e i c (R) is weaker 
than U g (R) oc R 6 in Eq.(l), as E el (R) ~ WR 2 for R « K/W, and E el (R) ~ KR for 
R » K/W [22]. 

We turn now to the "pure bulk" elastic interactions. In the hybrid aligned film, n is 
determined by the balance of anchoring and elasticity. In the thick films (h » K/W, 
h » R), n at the LC-air interface is vertical. Tangential anchoring of n at the droplet 
surface leads to a point defect-boojum [22] close to the south pole, Fig.le and 2e. The 
distortions n (r) around each droplet can be approximated by an elastic dipole normal to 
the LC-air interface [13], p ~ (0,0, a z R 2 ), where a z is a constant. The bulk elastic dipole- 
dipole interaction is thus isotropic in the xy plane and repulsive, 

a 2 R A 

U elb , thick (r) = K^-. (3) 
The total pair potential, U = U e i c + U e ib t thick, has a clear minimum at 

C _„2_ D4r/\l/3 



e.g., rlf r = 13/xm for f d = l()- 10 N , a z = 0.2 and R = 3fim. For f d oc WR, Eq.(4) predicts 
r p«r = 7 ( aLCA R^K/W 2 ) 1/3 oc i? 2 / 3 , where 7 is a dimensionless constant. The trend is in 
agreement with the experiment, Fig.2d, despite the fact that Eq.(4) is derived from the 
pair potential. Generally, f e i might be a non-linear function of R and also might depend 
on r. A surfactant can be used to control the patterns, as it changes <7 LC a, anchoring 
at interfaces, and thus j e \. Sometimes, one observes the patches of a hexagonal "dense" 
pattern with r eg:d (2.2 — 2.7)/? that coexist with the regular hexagonal patterns of a 
larger period r eq (3 — 4) R (but the same R). These dense patches are not understood 
yet, and might be provoked by ionic and flexoelectric charge, by the r-dependence of f e i, 
surfactant contamination, etc. 

In the thin films, n becomes progressively horizontal at the LC-air interface as the 
anchoring coefficient W air there is smaller than W [26]. When h decreases to 2R and K/W air , 
the boojum shifts to the side, Fig.3b, giving rise to the in-plane elastic dipole component, 
Px = o: x R 2 , and an anisotropic interaction contribution 

U elb , thin (r, 9) = K^- (1 - 3 cos 2 9) , (5) 

that depends on the angle 9 between the radius-vector connecting the two droplets and 
the direction of p x 's. The chain formation can thus be explained similarly to the chaining 
of droplet-hedgehog pairs entirely in the LC bulk [13]: U e ib t thin is minimum when 9 = 0. 
The analogy is not complete, as our problem is of a lower symmetry and requires one to 
consider h-, W a i r -, and W^-dependences of n(r), etc. For example, if the film is so thin that 
n is parallel to the LC-air interface [26], then n(r) around the droplets might become of a 
quadrupolar type with two boojums along the horizontal axis. 

To conclude, the particles (glycerol droplets) at the LC-air interface are capable of both 
repulsive and attractive interactions that lead to ordered patterns of hexagonal and chain 
types. The director distortions in LC give rise to a vertical elastic force balanced by surface 
tension; the deviation of the LC-air interface from the horizontal plane results in attraction 
of particles. The phenomenon is general, as we observe hexagonal patterns not only for the 

6 



liquid droplets but also for micron-size polymer spheres at the nematic-air interface; these 
results will be published elsewhere. The elastic-capillary mechanism can also shed some 
light at the observations of hexagonal patterns claimed to form in the bulk of the hybrid 
aligned films [18-20]. If the particles were indeed in the bulk, where n is tilted, then the 
elastic interaction potential would be anisotropic and thus generally inconsistent with the 
hexagonal order. The elastic model [18] does not consider the angular dependence of the 
potential for tilted n and does not explain the hexagonal ordering. A pure elastic model 
would not be satisfactory even for thick films, in which n around the droplets is vertical, 
as in this case the interaction would be only repulsive, Eq.(3). Note that the experiments 
[18-20] relied on the standard (non-confocal) microscopy with a poor z-axis resolution. The 
particles might have been located at the interface rather than in the LC bulk, which would 
bring the observations [18-20] into the category of elastic-capillary phenomenon proposed 
in this work. 

We acknowledge helpful discussion with M. Kleman, V. Pergamenshchik, S.V. 
Shiyanovskii, P. Tomchuk, and H. Yokoyama. The work was supported by National Science 
Foundation Grant DMR-0315523 and STCU grant #2025. ODL thanks the participants 
of the "Geometry and Materials Sciences" Program at the Aspen Center for Physics for 
discussions. 
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Figure captions. 

FIG.l. Two-channel FCPM textures of the vertical cross-section of the LC film: (a,b) 
a single glycerol droplet at the LC-air interface; (c) a raw of droplets that is a part of a 
hexagonal pattern at the LC-air interface; (d) scheme of forces and interfacial tension vectors 
at the LC-air interface with a strongly exaggerated meniscus slope; (e) 3D director field in 
the LC layer. The LC layer (glycerol droplet) is manifested by the high (low) intensity of 
fluorescent light in parts (a,c); the contrast is opposite in (b). 

FIG. 2. Optical microscopy pictures of hexagonal structures at the free surface of the 
thick nematic layer, h ~ 60/im, formed by droplets of different average diameter: (a) 2R m 
7/j,m, (b) 2R l/im. (c); arrangement of droplets after the surfactant is added to the 
nematic layer; (d) droplets separation r eq vs. 2R; the data are compared to the dependence 
r pair = 1 i K a LCA R 2 K/W 2 ) 1 ^, where 7 = 0.43, K = lO" 11 ^, a LCA = 3.8 x 10" 2 J/m 2 , 
W = 1CT 5 J/m 2 , see the text; (e) schematic drawing of the director field in the top part of 
the LC layer. 

FIG. 3. Optical microscopy textures of glycerol droplets forming chains at the free surface 
of a thin LC layer, h (7 — 10) jim (a) and the corresponding director field scheme (b). 
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